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A method to calculate a corrected effective diameter (Di) to more accurately 
estimate the dose received by a patient in chest computed tomography (CT) 
examination had been previously proposed. However, the discrepancy 
between Dmm and water-equivalent diameter (Dw) is still relatively high (i.e. 
about 6%). Furthermore, the method is still performed manually, so it is 
laborious and time-consuming. This study aims to improve the corrected 
effective diameter with bone correction (Dert®°") and to automatically calculate 
it. The automated Deri*°" was calculated as the square root of the product of 
these corrected AP and LAT diameters. The approach was implemented on 30 
patients who had undergone chest CT examination with the standard protocol. 
The results show that the correlation between the Deti°°" and Dw is R7=0.93 
with no statistical difference (p>0.05). The automated Derr" is 3.1% lower 
than Dw. While the Dum is 10.5% higher than Dw and both are statistically 


IndoseCT different (p<0.05). In conclusion, the new Dert*" was introduced and the result 

Size-specific dose estimate obtained was closer to Dw than Dm. This method is simple enough to be used 

Water-equivalent diameter as an alternative method to accurately estimate Dw for radiation dose 
estimation in clinical chest CT scanning. 
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1. INTRODUCTION 

Computed tomography (CT) is one of the medical imaging modalities widely for early diagnosis of 
disease [1], [2]. Although it has enormous benefits, CT delivers relatively high radiation dose to the patient [3], 
[4]. It was reported that the dose received by patients in CT examination is about 15 mSv in an adult and 30 
mSv in neonates for every single CT examination [5], whereas examinations with mammography, conventional 
X-ray, and nuclear medicine deliver radiation doses of 0.3-0.6 mSv, 5-10 mSv, and 2-5 mSv, respectively [6]. 
Moreover, the dose delivered by a single CT examination is greater than the annual radiation dose received 
from natural radiation sources such as radon and cosmic radiations (1-10 mSv) [7]. 

The United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) reported 
that high radiation doses pose a higher risk of cancer in the general population [8]. Consequently, the topic of 
CT dosimetry is always of high interest [9], [10]. Radiation dose monitoring in CT examination is still based 
on the volumetric computed tomography dose index (CTDI,o1) and the dose-length product (DLP). However, 
these dose indices do not consider the size of patients [11]. It was reported that for the same CTDI,o1 value, the 
patient dose increases with decreasing patient size [12], [13]. Therefore, although the CTDI,.1 and DLP indices 
are accurate in estimating the radiation output of a CT scanner, they are less accurate in assessing the dose 
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received by the patient [14], [15]. The American Association of Physicists in Medicine (AAPM) in 2011 
introduced another dose index that considers the patient size, the size-specific dose estimates (SSDE) [16]. The 
size of a patient is initially represented by the effective diameter (Derr) [15]. However, it is known that Der does 
not consider the different attenuation properties of the various body parts. Differences in attenuation affect the 
absorption of radiation which impacts the dose received by the patient [17]. 

Water-equivalent diameter (D,,) was subsequently introduced by AAPM in task group 220 to 
incorporate both the size and attenuation properties of a patient’s body [18]. Automated methods to obtain the 
value of D, have been proposed by many researchers recently, such as Anam et al. [19], Ozoyka et al. [20], 
Gharbi et al. [21], and Juszczyk et al. [22]. Calculation of Dy needs specific software since it involves 
calculation of the average CT number and the area of the patient’s image [23]. Not all hospitals have such 
software so an alternative method for estimating Dy is needed. 

Mihailidis et al. [24] introduced a method to calculate the value of corrected effective diameter as an 
alternative to calculate the Dy in thoracic CT examinations. The lateral (LAT) thickness is corrected by the 
average relative electron density of the lung (pe=0.3) to overcome the presence of the lung, so that the results 
obtained are close to the value of Dy. However, the values obtained from this method are significantly different 
(i.e. about 6%) from the value of Dy because the presence of bone is not considered. Furthermore, this method 
is carried out manually using an electronic ruler and involves many steps. The purpose of this study is to 
improve the effective diameter correction with both lung and bone to obtain the result closer to the Dy and to 
propose an automated method for calculating it. 


2. METHOD 
2.1. The diameters calculation 

If we assume that the cross section of the patients is circular or elliptical, we can estimate the effective 
diameter of patients from the magnitude of the diameters in the anterior-posterior (AP) and lateral (LAT) 
directions (1). 


Depp = VAP X LAT (1) 


But, in fact, the patient’s cross section is neither elliptical or circular [25]. In addition, Derr does not consider 
the attenuation properties of body. 

The AAPM proposed the water-equivalent diameter (Dy) to incorporate both the size and attenuation 
properties of the patient’s body [18], as seen in (2). 


Dy =2x | (Oman +1) x [= (2) 
Where ROInmean is the mean CT number of the region of interest and Aror is total area of ROI of the axial image. 
In the Dy calculation, specific software is needed to calculate the average CT number and the area of the 
patient’s image. 

Mihailidis et al. [24] proposed a different approach for considering the lung tissue in the image. They 
corrected the lateral diameter with the relative electron density of lung tissue (p.=0.3), naming it the LAT ef", 
as seen in (3). It was hoped that this value would be equivalent to Dy [24]. 


Dun = [AP. LATOR (3) 


We postulate that to achieve more accurate results, the correction should not only be made in the LAT direction 
but also in the AP direction. Correction must account for the presence of bone as well as the presence of lung 
tissue. Therefore, we propose the corrected effective diameter (Defi°°") as seen in (4). 


Ce APL AT Core (4) 


In addition, we propose an algorithm to automate the calculation of Deg" to make it easier to 
implement in a clinical setting. The steps are shown in Figure 1, while the example from patient image is shown 
in Figure 2. First, we open the original image as seen in Figure 2(a). Then, we segmented the patient border 
automatically and converted the CT values outside the patient to a pixel value of 0 Figure 2(b). Then we 
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detected the presence of lung and bone. The lung was detected with a threshold of -500 HU and bone with a 
threshold of +100 HU. CT values lower than -500 HU were converted to a pixel value of 1. CT values in the 
range -500 HU<x<+100 HU were converted to a pixel value of 2. All other CT values were converted to a 
pixel value of 3. We then determined the center position of the patient and determined the diameters in AP and 
LAT from the central image Figure 2(c). In the lung tissue (i.e., with pixel value of 1), the AP and LAT 
diameters were corrected with the average relative electron density of lung (pe=0.3), and in the bone (i.e., pixel 
value of 3) corrected with the average relative electron density of bone (p-=1.8). Finally, the Der °°" was 
calculated using (4) and as seen in Figure 2(d). The algorithm was effectively integrated into a graphical user 
interface (GUI) of IndoseCT 20b as depicted in Figure 3. 
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Figure 1. Flowchart for automatic corrected effective diameter (Der®") calculation 


(c) (d) 


Figure 2. The segmentation process of Des®", (a) original image of patient, (b) the image after binarization, 
(c) the central position of LAT and AP diameters to estimate the central effective diameter (Derr), and 
(d) Result of patient image with changes in lung pixel value to 1, bone to 3, and otherwise to 2 
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Figure 3. The screenshots of IndoseCT 20b for an automatic calculation of corrected effective diameter in 
thoracic region 


2.2. Patient images 

In this study, we evaluated 30 patients who undergone chest CT examination. The examinations were 
performed using a 128-slice multi-detector CT scanner, the Toshiba Aquilion 128. The patients were scanned 
with standard imaging protocol using a voltage of 120 kVp, 3D tube current modulation (TCM), a pitch of 
1,438 and a collimation beam width of 640.5 mm in Ken Saras Hospital Semarang, Indonesia. 


3. RESULTS AND DISCUSSION 

The results of the calculated diameters (Dm, De®", and Dy) are shown in Figure 4. The average and 
standard deviation of the diameters (Defi®"", Dut, and Dy) are listed in Table 1. The Der" is slightly smaller 
than Dy. The percentage differences between Dm, Defi®°", and Dy are shown in Table 2. The differences 
between De" and Dy is 3.1 % and they are not statistically significant (p>0.05). However, Dum is 10.5 % 
different from D,,. The Dum is statistically different (p<0.05) from Dy. 

The relationship between the diameters is shown in Figure 5, with the automated Dy and Deg°" in 
Figure 5(a), and the automated D,, and the Dm in Figure 5(b). The automated De" is linearly correlated 
(R7=0.93) with Dy. The correlation between Dy and Dum is R7=0.95. 

The purpose of this study was to improve the corrected effective diameter with bone correction 
(Defi®") and to make it is automatically carried out. Previously, Mihailidis et al. [24] proposed a manually 
corrected Des (Dm) and the results were expected to close to Dy. However, their results were still 6.0 different 
from Dy [24] because the correction was only performed in the LAT direction and for lung only. Furthermore, 
their method is laborious and time-consuming. We have proposed corrections in two directions (AP and LAT) 
and inclusion of bone in addition to lung. We subsequently develop a software to automatically calculate the 
corrected Derg (Deri®°") which works quickly and. The corrections are crucial to obtain a value of Des®" close to 
D, so that the calculation of SSDE is more accurate. 

In Table 2, the value of automated Der" is only slightly lower (3.1%) than the value of Dy. This is 
because the calculation of Deg" considers lung and bone of the thorax. The presence of lung and bone in 
thorax are corrected with the average relative electron density of lung (pe lung=0.3) and bone (pe bone=1.8) in 
both the LAT and AP directions. In this study, the Dymis 10.48% lower than Dy. The larger difference between 
Dm and Dy is because Dm corrected the image with the attenuation of the lung in only one direction (only 
in the LAT diameter), while the value of Dy was corrected in two directions [24]. Ignoring bone in AP 
dimension in the Dv method made the values obtained far lower than Dy. 

The main finding of our study is that the automated corrected effective diameter can be close to water- 
equivalent diameter if we use a correction for the AP diameter as well as considering the presence of bone. The 
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automated De" can be used as a viable alternative since it has the close value to Dy. However, automatic 
detection of lung and bone can be tricky because it depends on the pixel thresholds used. In this study, we used 
-500 HU to detect the lung and +100 HU to detect the bone. The thresholds maybe different for other images. 
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Figure 4. Results of De", Dum, and Dy in thoracic examinations 


Table 1. The Dum, Dert®°", and Dy values 


Parameter Doo De" Dy 
Mean (cm) 20.82 22.93 23.21 
Std deviation (cm) 3.22 3.37 3.21 
Max (cm) 27.39 30.78 30.1 
Min (cm) 15.39 17.17 16.78 
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Table 2. Percentage differences the Der" and Du from Dw 
Parameter Mean Std Deviation Max Min _ p value 


Der" (%) 3.06 3.61 9.21 0.09 0.72 
Dou (%) 10.48 2.22 16.76 2.42 <0.01 
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Figure 5. The relationships between (a) Dum and Dy, and (b) De®" and Dy, in thoracic examinations 
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4. CONCLUSION 

The new Deg" was introduced and the result obtained was closer to Dy than Dum. A software to 
automatically estimate Dep" has been developed and tested. It can be operated easily and quickly. The 
corrected effective diameter gives the close result to Dy (i.e. the difference is about 3%) and it can be an 
alternative to measuring Dy in daily examinations in routine clinical chest CT. 


ACKNOWLEDGEMENTS 
This work was funded by the World Class Research University (WCRU), Diponegoro University, No. 
118-08/UN7.6.1/PP/2021. 


REFERENCES 

[1] J. Lee et al., “Neural network-based method for diagnosis and severity assessment of Graves’ orbitopathy using orbital computed 
tomography,” Scientific Reports, vol. 12, no. 1, Jul. 2022, doi: 10.1038/s41598-022-16217-z. 

[2] A. Ismail et al., “Clinical and chest computed tomography features of patients suffering from mild and severe COVID-19 at Fayoum 
University Hospital in Egypt,” PLoS ONE, vol. 17, no. 7, Jul. 2022, doi: 10.1371/journal.pone.0271271. 

[3] R. Huang, X. Liu, L. He, and P. K. Zhou, “Radiation Exposure associated with computed tomography in childhood and the 
subsequent risk of cancer: A meta-analysis of cohort studies,’ Dose-Response, vol. 18, no. 2, Apr. 2020, doi: 
10.1177/1559325820923828. 

[4] C. Anam, F. Haryanto, R. Widita, I. Arif, and G. Dougherty, “The evaluation of the effective diameter (Deff) calculation and its 
impact on the size-specific dose estimate (SSDE),” Atom Indonesia, vol. 43, no. 1, pp. 55-60, Apr. 2017, doi: 10.17146/aij.2017.617. 

[5] D. J. Brenner and E. J. Hall, “Computed tomography — an increasing source of radiation exposure,” New England Journal of 
Medicine, vol. 357, no. 22, pp. 2277-2284, Nov. 2007, doi: 10.1056/nejmra072149. 

[6] L. Yu etal., “Radiation dose reduction in computed tomography: techniques and future perspective,” Imaging in Medicine, vol. 1, 
no. 1, pp. 65-84, Oct. 2009, doi: 10.2217/iim.09.5. 

[7] C. H. McCollough, L. Guimaraes, and J. G. Fletcher, “In defense of body CT,” American Journal of Roentgenology, vol. 193, no. 
1, pp. 28-39, Jul. 2009, doi: 10.2214/AJR.09.2754. 

[8] | United Nations Scientific Committee Effects of Atomic Radiation, “Report to the general assembly scientific Annex A:Volume I,” 
USA: United Nations, 2014. [Online]. Available: https://www.unscear.org/unscear/en/publications/2020_2021_1.html. 

[9] S. P. Raman, M. Mahesh, R. V. Blasko, and E. K. Fishman, “CT scan parameters and radiation dose: Practical advice for 
radiologists,’ Journal of the American College of Radiology, vol. 10, no. 11, pp. 840-846, Nov. 2013, doi: 
10.1016/j.jacr.2013.05.032. 

[10] S. Maxwell et al., “How have advances in CT dosimetry software impacted estimates of CT radiation dose and cancer incidence? 
A comparison of CT dosimetry software: Implications for past and future research,” PLoS ONE, vol. 14, no. 8, Aug. 2019, doi: 
10.137 1/journal.pone.0217816. 

[11] C.H. McCollough, S. Leng, L. Yu, D. D. Cody, J. M. Boone, and M. F. McNitt-Gray, “CT dose index and patient dose : They are 
not the same thing,” Radiology, vol. 259, no. 2, pp. 311-316, May 2011, doi: 10.1148/radiol.11101800. 

[12] W. Huda and F. A. Mettler, “Volume CT dose index and dose-length product displayed during CT: What good are they?,” 
Radiology, vol. 258, no. 1, pp. 236-242, Jan. 2010, doi: 10.1148/radiol.10100297. 

[13] J. A. Christner, N. N. Braun, M. C. Jacobsen, R. E. Carter, J. M. Kofler, and C. H. McCollough, “Size-specific dose estimates for 
adult patients at CT of the torso,” Radiology, vol. 265, no. 3, pp. 841-847, Dec. 2012, doi: 10.1148/radiol. 12112365. 

[14] C. Anam, F. Haryanto, R. Widita, I. Arif, and G. Dougherty, “A fully automated calculation of size-specific dose estimates (SSDE) 
in thoracic and head CT examinations,” Journal of Physics: Conference Series, vol. 694, no. 1, Mar. 2016, doi: 10.1088/1742- 
6596/694/1/012030. 

[15] “The measurement, reporting, and management of radiation dose in CT,” American Association of Physicist in Medicine, Jan. 2008. 
doi: 10.37206/97. 

[16] “Size-specific dose estimates (SSDE) in pediatric and adult body CT examinations, TG-204,” American Association of Physicist in 
Medicine, 2011. 

[17] J. Wang, X. Duan, J. A. Christner, S. Leng, L. Yu, and C. H. McCollough, “Attenuation-based estimation of patient size for the 
purpose of size specific dose estimation in CT. Part I. Development and validation of methods using the CT image,” Medical 
Physics, vol. 39, no. 11, pp. 6764-6771, Oct. 2012, doi: 10.1118/1.4754303. 

[18] C. McCollough et al., “Use of water equivalent diameter for calculating patient size and size-specific dose estimates (SSDE) in 
CT,” American Association of Physicist in Medicine, AAPM report, pp. 6-23, 2014. 

[19] C. Anam, F. Haryanto, R. Widita, I. Arif, and G. Dougherty, “Automated calculation of water-equivalent diameter (DW) based on 
AAPM task group 220,” Journal of Applied Clinical Medical Physics, vol. 17, no. 4, pp. 320-333, Jul. 2016, doi: 
10.1120/jacmp.v17i4.6171. 

[20] £. Ozsoykal, A. Yurt, and K. Akgiingdr, “Size-specific dose estimates in chest, abdomen, and pelvis CT examinations of pediatric 
patients,” Diagnostic and Interventional Radiology, vol. 24, no. 4, pp. 243-248, Jul. 2018, doi: 10.5152/dir.2018.17450. 

[21] S. Gharbi, S. Labidi, and M. Mars, “Automatic brain dose estimation in computed tomography using patient dicom images,” 
Radiation Protection Dosimetry, vol. 188, no. 4, pp. 536-542, Jul. 2020, doi: 10.1093/rpd/ncaa006. 

[22] J. Juszczyk et al., “Automated size-specific dose estimates using deep learning image processing,” Medical Image Analysis, vol. 
68, Feb. 2021, doi: 10.1016/j.media.2020.101898. 

[23] C. Anam et al., “An improved method for automated calculation of the water-equivalent diameter for estimating size-specific dose 
in CT,” Journal of Applied Clinical Medical Physics, vol. 22, no. 9, pp. 313-323, Sep. 2021, doi: 10.1002/acm2. 13367. 

[24] D. Mihailidis, V. Tsapaki, and P. Tomara, “A simple manual method to estimate water-equivalent diameter for calculating size- 
specific dose estimate in chest computed tomography,” British Journal of Radiology, vol. 94, no. 1117, Jan. 2021, doi: 
10.1259/bjr.20200473. 


Indonesian J Elec Eng & Comp Sci, Vol. 31, No. 1, July 2023: 189-196 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 ) 195 


[25] C. Anam, F. Haryanto, R. Widita, and I. Arif, “Automated estimation of patient’s size from 3D image of patient for size specific 
dose estimates (SSDE),” Advanced Science, Engineering and Medicine, vol. 7, no. 10, pp. 892-896, 2015, doi: 
10.1166/asem.2015.1780. 


BIOGRAPHIES OF AUTHORS 


Dr. Choirul Anam © &:4 3 & completed his Ph.D. in Physics Department, Bandung 
Institute of Technology (ITB). He received Master degree from University of Indonesia 
(UD) and the B.Sc. degree from Diponegoro University (UNDIP). He is currently working 
as a Lecturer and Researcher at the UNDIP. His research interests are medical image 
processing and dosimetry for diagnostic radiology, particularly in CT. He has authored and 
co-authored over 150 papers. One of his papers published in the Journal of Applied Clinical 
Medical Physics (JACMP) had been awarded as the Best Medical Imaging Physics Article 
in 2016. He was also recognized as an Outstanding Reviewer for the Physics in Medicine 
and Biology (PMB) in 2018 and for the Biomedical Physics and Engineering Express 
(BPEX) in 2019. He is developer of two software, i.e. IndoseCT (for calculating and 
managing radiation dose of CT) and IndoQCT (for measuring CT image quality). He can 
be contacted at email: anam @fisika.fsm.undip.ac.id. 


Riska Amilia © 4 1D received her B.Sc. in Diponegoro University. Her research 
interests are radiation physics, radiation protection, CT dosimetry, and medical imaging. 
She can be contacted at email: riskamilia0 @ gmail.com. 


Prof. Wahyu Setia Budi © (S he is a professor in Department of 
Physics, Diponegoro University. His research of interests are radiation physics, nuclear 
physics, optics, and X-ray imaging. He can be _ contacted at email: 
wahyu.sb @fisika.fsm.undip.ac.id. 


Prof. Heri Sutanto ‘ Ef (D is a lecturer in materials physics at Diponegoro 
University. His expertise is thin films, nanomaterials, and medical materials. He can be 
contacted at email: herisutanto @live.undip.ac.id. 


An improvement of the computational effective diameter measurement in thoracic ... (Choirul Anam) 


196 0 ISSN: 2502-4752 


Zaenul Muhlisin © £4 B " is a lecturer in medical physics at Diponegoro University. 
His expertise is X-ray imaging and radiation protection. He can be contacted at email: 
muhlisin @lecturer.undip.ac.id. 


Ariij Naufal © ki 1S is a graduate student of the physics master's program at 
Diponegoro University. He is the developer of IndoQCT. He can be contacted at email: 
ariij.2019 @fisika.fsm.undip.ac.id. 


Prof. Geoff Dougherty ONES graduated with First Class in Physics at Manchester 
University, did a post-baccalaureate teaching certificate at Leeds University, taught at a 
community college and then completed a Ph.D. on DNA-drug interactions at Keele 
University. This was followed by a post-doc period investigating closed-circular DNA at 
the Swiss Federal Institute of Technology, Zurich. He taught at the Science University of 
Malaysia; did research on DNA-drug binding using electron spin resonance at Monash 
University, Australia; and taught at the University of the South Pacific, Fiji and Oxford 
Brookes University. After that he moved to Kuwait as Professor of Radiologic Sciences 
(1992-2002). He moved to California in 2002, where he is currently Professor of Applied 
Physics and Medical Imaging. He applies image analysis and pattern recognition 
techniques to medical images from a variety of modalities in an effort to extract the 
maximum quantifiable diagnostic and prognostic information from them. His three 
Fulbright experiences (two as Senior Scholar and one as a Specialist) have resulted in many 
new collaborations. He can be contacted at email: Geoff.Dougherty @csuci.edu. 


Indonesian J Elec Eng & Comp Sci, Vol. 31, No. 1, July 2023: 189-196 


